clear
clc
N=40
maxscore=0
length_days=90
T=360; %number of periods
Svalues=unique([unique(maxscore);max(maxscore)+[0.002:0.002:0.05]']);
nScore = size(Svalues,1);
lambda1=0.75
lambda2=1-lambda1;
sigmaS=0.025;
sigmaT=0.025;
nSim=500;
mergers_data=5;

q_sp=exp(-3*Svalues-1)./(1+exp(-3*Svalues-1));
q_team=exp(-3*Svalues-0.5)./(1+exp(-3*Svalues-0.5));